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Using  a  general  method  for  deriving  identities  for  random  variables,  i 
find  a  number  of  new  results  involving  characteristic  functions  and 
generating  functions.  The  method  is  simply  to  promote  a  parameter  in 
an  integral  relation  to  the  status  of  a  random  variable  and  then  take 
expected  values  of  both  sides  of  the  equation.  Results  include  formulas 
for  calculating  the  characteristic  functions  for  x2,  /x,  l/x,  x2  +  x,  R2  = 
x2  -  g2,  etc.  ip  te^ms  of  integral  transforms  of  the  characteristic 
functions  for  x  and  (x.  u),  etc.  Generalizations  to  higher  dimensions  can 
oe  obtained  us  mg  the  same  method.  Expressions  for  inverse/fractionai 
moments.  E(n?}.  etc.  are  also  presented,  demonstrating  the  method. 


. .  , .  i  ,  rtUu’..  -w  .  L’  i  h 

as  s  well  known,  it  is  sometimes  easier  to  study  a  process  using 
‘cans  forms  of  the  relevant  probability  distributions.  Such  transforms 
•r c  ude:  the  characteristic  function,  C(o>),  and  the  moment  generating 
•'unction,  M(G).  for  general  random  variables;  the  probability  generating 
•  jnrnon.  G(z;,  ’or  integer  valued  random  variables,  and  the  Laplace 
’'■jrsUrr:'  of  the  probability  density  function,  $As).  for  non-negative 
valued  random  variables.  They  often  allow  one  to  I)  simplify 
manipulations  involving  convolutions  of  probability  distributions  arising 
from  consideration  of  sums  of  random  variables  and  more  complicated 


compound  and  branching  processes,  and  2)  apply  powerful  methods  from 
complex  analysis  and  integral  transform  theory  to  the  solution  of 
differential-difference  equations  which  arise  in  the  study  of  probability 
and  stochastic  processes,  and  in  the  analysis  of  the  analytic  behavior  of 
those  solutions.  The  value  of  techniques  for  manipulating  such 
transforms  and  of  "methods  for  constructing  new  characteristic 
functions  out  of  given  ones"1  is  well  known.  In  fact,  the  theory  of 
probability  "depends  to  a  large  extent  on  the  method  of  characteristic 
functions"2.  The  methods  presented  here  may  further  aid  in  the 
interpretation  of  complicated  characteristic  functions  and  facilitate  the 
identification  of  independent  processes  which  contribute  to  the  result ' 
'.see  e.g.  reference  [3]).  Apart  from  their  usefulness  in  probabilistic 
applications,  our  results  also  provide  another  means  of  generating  new 
integral  identities  from  old  ones. 

By  promoting  a  parameter  in  an  integral  expression  to  the  status  of  a 
random  variable  (r.v.)  and  then  taking  expected  values  of  both  s'des  of 
the  equation,  a  number  of  interesting  relations  involving  characteristic 
functions,  generating  functions,  etc.  are  found.  In  general,  while  there 
is  no  guarantee  that  the  resulting  integrals  can  be  evaluated  in  ciused 
form  for  a  particular  distribution  of  interest,  the  expression  may  be 
Helpful  in  numerical  work.  An  analogous  technique  for  generating 
'dennties  nvolvmq  operators  in  Hilbert  space,  matrices,  etc.  has  been 
useful ;.y  employed  in  physics  (e  g.  see  the  Appendix,  below),  in  tne 
probability  context  sirmlar  methods  have  long  been  used  to  solve 
problems  by  averaging  conditional  results  over  me  conditioning  variable. 

2 


»tn. 


A  number  of  identities  are  presented  in  Sections  ii  -  V  demonstrating 
me  method  of  derivation.  Some  examples  of  calculations  using  these 
^entities  are  then  carried  out  in  Section  V!.  Finally,  in  Section  vn  we 
comment  on  the  generality  of  the  method. 


!!.  Relations  involving  the  Souare  of  R,  V,  s 

a^  Consider  the  well-known  integral  expressing  the  normalization  of  a 
Normal  (aka  Gaussian)  distribution,  in  which  x  is  an  arbitrary  constant, 

oo 

!'exp[  -(£  -  x)2/(2o2)  1  d£  =  1  . 

*°°  /(2m  a2) 

Change  variables  according  to.  £  — >  £/(22f).  l/(2a2)  -->  i2f,  an 


CC 

1  e > d f - 1 c -J 7 1 4 )  *  i>Cx  i  d£  =  exp[  itfx^  ]  . 
*  /(-i4tt20 


(1 


Now  promnte  x  to  be  a  real  random  variable  and  take  expected  values  of 
Doth  sides  of  the  equation,  assuming  that  the  implicit  interchange  of 
:rliirs  :f  integration  is  justified,  i.e.  that  [  and  E{  }  commute, 

00 

=  CX?(V)  -  (2) 


/[  1  Jexp[-i£2/(4tf)l  Cx(£) 


oo 


!  I expf~i£2/(4ff)  *  ifx  -i62/(4S)  lifyl  dfde  =  exp[  i3x2  ♦  iSy2  ]  . 

[  -!4tt/(^8)  1 

Aga:r:  consider  x,  y  to  be  r.v.'s  and  take  expected  values  of  both  sides, 

00 

i/lW(3S)  1  J JexpH£2/(4ff)  -if2/(4S)]  Cx  y(£.e)  d£de  (3) 

-oo 

CX2,y2(2f,S)  . 

If  we  now  let  8  =  Z  we  have, 

00 

i/l«ir*l  J/expHtt2  *  «J)/(4*)1  Cx _„«.()  d?dt  (4) 

-00 

=  Cr2(*)  . 

where-  R2  =  X2  *  V2.  This  can  be  generalized  further  to  3  or  more  r.v.’s 
■n  an  analogous  manner, 

_  Multiply  Eg.  (1)  by  exp(  itfx)  and  take  expected  values  to  obtain  the 
characteristic  function  for  X2  ♦  x. 

00 

A  i/(4*a)  1  Jexpl-i£2/(4ff)l  Cx(£  ♦  #)  d£  =  Cx2*x(#)  .  (5) 

-oo 

Again,  it  is  clear  that  this  may  be  generalized  further. 


a_  Consider  the  definite  integral  (e.g.  reference  [6]  p.  34 1). 

00 

f  exp  [  -a/£2  -b^2  1  d£  =  /[tt /(4b)]  exp(-2/(ab)l  .  (6) 

0 

Let  a  -->  b  -->  s2/4  to  obtain  the  identity, 

00 

\  exp  [  -x/£2  -s2£2/4  i  d£  =  /(tt)  /s  exp[-s/ x]  .  (7) 

o 

Now,  promote  x  to  be  a  non-negative  r.v.  and  average  over  x,  to  obtain 
the  Laplace  transform  of  the  pdf  of  /x, 

00 

s//ir  J  exp  [-s2£2/4  1  2,x(l/*;2)  =  2, /x(s)  .  (8) 

o 

Alternatively,  a  similar  integral  on  p.  399  of  reference  [6]  allows  one  to 
express  SL yx(s)  in  terms  of  the  characteristic  function.  Cx  . 

±  To  obtain  the  Laplace  transform  of  the  pdf  for  the  r.v.  i/x,  i.e. 

2  1  '  y  ,  given  SLX  ,  consider  the  integral 10, 

X) 

i  pxd[  -a£  i  J0(b/£)  dc,  =  l/a  exp[  -  b2/(4a)  1. 

!J 

Mult.piu  Doth  sides  by  a  and  change  the  parameters  b  --s  2v's,  a  v 
the  i.atter  a  non-negative  r.v.,  to  obtain. 


00 


{  Ef  x  exp[  -xt  } }  •  J0[  2/(sC)  ]  d£  =  E{  expf  -  s/x  !  }. 
o 

In  terms  of  the  Laplace  transform,  this  is, 

00 

(-)  J  -  Jol  2/(sO  ]  dS  =  at,/x(s).  (9) 

0 

iv.  Identities  for  Non-Standard  Moments  and  Averages 
A.  Consider  the  elementary  integral,  where  x  is  just  a  parameter, 

00 

j  sn‘ 1  exp  [  -x  s]  ds  =  (n- 1 )!  /  xn  . 
o 

Now  promote  x  to  be  a  non-negative  r.v.,  whose  pdf  falls  off  sufficiently 
rap.diy  as  x  -->  C.  (e.g.  an  Erlang(n+1)  )  and  take  expected  values  t  x. 


oo 

1  /(n- 1 )!  /  s"'1  Xx(s)  ds  =  E{l/xn  }  (10a) 

o 

ettmq  x  -->  (x  +  A)  leads  immediately  to  the  identity, 


CO 


An- i)i  i  sn  1  exp(-sA')  SL. (s)  ds  =  E{l/fx+A)p 


(10b) 


idenf'fes  for  the  Laplace  transform  could  also  be  written  in  terms  of 


<*A 


the  moment  generating  function,  when  it  exists.  Analogous  results  for 
the  moment  generating  function  were  also  derived  in  references  [4],  and 
25!  using  methods  similar  to  the  above.  Those  references  also  contain 
additional  applications  of  this  result. 

EL  Consider,  now,  the  integral. 


2  |  exp  [  -a  t2]  dt  =  /(tt/a) 


let  a  — >  x,  a  non-negative  r.v.,  and  take  expected  values, 


1 1 


00 


00 


2//*  J  Cx(it2)  dt  -  2//tt  /ix(t2)  dt  -  E{l//x}  (12) 


Taking  the  change  of  variable  to  y  =  t2  in  Eq.(12)  results  in, 


1/2-1 


&x(y)  dy  =  E(’//x) 


'-■is  can  be  recognized  as  a  fractional  integration  of  order  1/2  of  the 
ap'.ace  transform  (or  nGF).  Some  of  the  other  moments  in  this  section 
:ar.  also  be  written  as  fractional  integro-differentiations  of  moment 
;*ner.3t.r.g  :  .notions  or  Laplace  transforms.  This  fact,  as  well  as  other 
sions  -anc  re.ated  references)  are  discussed  in  references  [7  -  9]. 


-  .  *  - 


•■**'  a..,.  aV  ■*"*•*-*  -UnTi*  rffrhti  r^>ihv 


C  Multiply  Eq.  (M)  by  a,  let.  a  *  and  average,  to  obtain, 

>0 

Z/'/tt  {  E(x  exp(-xt2]}  dt  =  El/x). 
o 

or,  switching  to  moment  generating  functions  instead  of  Laplace 
transforms  for  this  result  (either  could  be  used  here), 

00 

2//it  /  rrx(-t2)  dt  =  e{/x}.  (13) 

o 

This  can  be  generalized  to  obtain  a  formula  for  E(  xm+1/2  },  with  m  an 
integer,  in  a  straightforward  manner. 

CL  Consider  Lipschitz's  integral10  for  the  ordinary  Bessel  function  of 
zeroth  order,  J0, 

•j*0 

f  exp[  -as  ]  JG(bs)  ds  -  l/v^a2  +  b2). 

Dromote  a  -->  x,  a  non-negative  random  variable  and  take  expected 
values  cr  both  sides  to  obtain, 

00 

[  £x(s)  J0(bs)  ds  =  E(  l//(x2  ♦  b2)  ).  (M) 

o 

Successively  differentiating  tn is  identity  wrt  the  parameter  b  produces 
a  family  of  similar  identities. 


L  Consider  one  form  of  Bessel's  integral  for  the  n^  order  ordinaru 
Bessel  function10, 

TT 

Jp(x)  -  1/(2tt)  f  exp[  -n  i  0  ♦  i  x  sine  ]  d6, 

-7T 

let  x  be  a  r.v.  and  average  over  ali  x,  to  obtain, 

7t 

E{  Jn(x)  }  =  i/(27t)  /  expl  -n  i  0  1  cx(sine)  de.  (15) 

-ir 

Clearly,  this  result  can  be  generalized  in  many  ways,  and  is  somewhat 
rem:n'=cent  of  the  well-known  formula. 

E'  Bfx)  }  -  i/(2tt)  |  H(u>)  Cx(o))  dm, 


■"here  mw''  is  tne  Fourier  transform  of  H(x).  The  latter  eouation  can,  m 
rhe  spirit  of  this  paper,  be  simply  derived  by  taking  expected  values  of  - 
:n  me  -ec  esentat'on  of  H(x)  as  the  Fourier  Transform  of  Hfmi. 

■  id-rtties  for  Probability  Generating  Functions 
A  Tons  der  again  the  well-known  integral  used  to  define  *he  gamma 


Tms  time  let  n  be  a  non-negative  mteger  valued  random  variable  and 


2veraqe  over  n. 


OO 


{  G(z)  expl  -sz]  dz  =  E{  n!  /  sn*  1  }. 


(16) 


in  part.cuiar.  when  s  =  1,  this  yields  E{  m  !.  when  it  exists,  i.e.  the 
:_3D,ace  transform  of  the  probability  generating  function,  evaluated  at 
-  ’  's  lust  t[ m  }  (cf.  the  factorial  moments  E{n(n- 1  V--(n-k- i )}  - 


).  For  non- integer  r.v.'s  we  can  obtain  a  corresponding 


expression  for  E{  f(x)  }. 


L  Consider  the  integral, 


Now  'et  n  oe  a  non-negative  integer  valued  r.v.  and  average. 

i 

J  G(z)  dz  -  E{  1/(n+ 1)  },  (17) 

o 

.vn";  a i so  follows  easily  from  the  power  series  definition  of  G(z;  and  is 
o.rect.y  analogous  to  the  usual  result  for  E{n}. 

^  ce*'  tne  integral  expressing  the  standard  result  for  the  even 
enta  :f  tne  Normal  distribution. 


CO 


(2n  -  1)11  cr2n  , 


?xpi  ~zi 


/(  2ttc2) 

where  the  double  factorial  symbol  means,  e.g.,  5!! 
averages  over  n  on  both  sides  of  the  equality. 


=  5-3- 1 ,  Again  take 


E{  (2n  -  1)!!  cf2n  }.  (18) 


■°°  S(2t«j2) 

and  when  o  =  1  we  have  E{  (2n  -  1)!!  }. 


Cl  Consider  the  two  integrals,  found  on  p.  369  of  reference  (61, 


:sm20]m  d9 


tt/2  •  (2m  -  1)h/  (2m)!!  , 


J  sine  [sin2e]m  d6  =  (2m)11  /  (2m  +  !)•'  . 

o 

.erring  m  be  a  r.v.  and  averaging  over  all  values  of  m  on  each  side  of  the 
above  equations,  we  obtain. 


2/it  \  G(  sin2e  )  de 


E{  (2m  -  I)!!  /  (2m)!!  }  . 


(19) 


r  *  '  'V 


E{  (2m)H  /  (2m  ♦  1)1!  }  . 


(20) 


/sine  G(  sin2e  )  de 

o 

respectively. 

Vi.  Some  Applications  of  the  Identities 
A.  If  x  has  a  Normal  distribution  with  zero  mean  then 
Cx(£)  =  expK2o2/2}. 

Putting  this  in  Eg.  (2)  and  performing  the  integration,  we  have, 

00 

Cx2(3)  =  /[  i/(47rt0  ]  |exp['i^2/(4^)  -  £2a2/2  ]  d£  . 

-00 

=  l//d  -  2i022f).  (21) 

This  is,  indeed,  recognized  as  the  characteristic  function  for  the  X2 
distribution  with  one  degree  of  freedom.  (Similarly,  if  x,  y  have 
independent  normal  distributions  with  the  same  value  of  the  variance, 
that  R2  =  x2  *  Y2  has  a  negative  exponential  distribution  follows 
tr’vially  from  Eq.  (4).  ) 

Now.  let  x  have  a  Normal  distribution  with  non-zero  mean,  p,  then 

C X(C) z  exp{  ip£  -  £2o2/ 2}. 

Substituting  this  in  Eq.  (2)  and  integrating,  yields. 


/[  i/(47T25')  J  Jexp(  ip£  -  i«C2/(42<)  -  l202/ 2  ]  d£  . 


Cx2(25')  = 

-00 

=  l//(1-  2io2#)  exp{  ip2S/(l  -  2i 022f)  },  (22) 

which  is  the  characteristic  function  of  an  offset  X2  distribution. 

tL  Calculate  E{l//x  }  where  x  has  an  exponential  distribution  with 
parameter  X.  In  this  case  the  characteristic  function  of  x  is, 

Cy(0  =  X/(X  -  i  £  ). 

‘-'ence.  substituting  this  m  Eq.  (12)  we  have. 

00  00 

2//n  f  Cx(it2)  dt  =  2//tt  J  X/(X  +  t2)  dt, 
o  o 

and,  using  a  standard  integral,  we  obtain, 

E(l//x  }  =  tt/X.  (23) 

mis  is  easily  verified  to  be  correct  by  3  direct  calculation.  E{  /x  }  is 
also  easily  verified  to  be  the  result  produced  by  Eq.  0  3). 

Insert  tne  Laplace  transform  of  the  pdf  for  an  exponential 
distribution,  X/(X  +  s),  into  Eq.  (14),  obtaining, 

OO 

Ei  I // (  x2  +  b2  )  )  -  J  X/(X  +  s)  ■  J0(bs)  ds  . 


T h i s  integral  is  tabulated  on  p.  685  in  reference  [6],  resulting  in, 

E(  l//(  x2  ♦  b2  )  }  =  X  7T/2  [  H0(bX)  -  N0(DX)  ].  (24) 

where  H0  and  N0  are  Struve  and  Neumann  functions,  respectively,  of 
zeroth  order  (N0  can  be  replaced  with  v0,  the  Bessel  function  of  the 
second  Kind).  For  example,  taking  b  -  4  and  X  =  1  (  H0(4)  =  .1350 land 
Yg(4)  =  -.01694  )s  we  find  for  the  exponential  distribution, 

E{  1//(  x2  >  D2  )  )  =  .2387  . 

which  is  easily  confirmed  by  direct  Gauss-Laguerre  integration  of  the 
left-hand-side. 

Cl  we  now  calculate  the  average  of  the  n^h  order  Bessel  function  when  x 
has  a  N(0,  c)  distribution  with  the  use  of  Eg.  (15).  After  inserting  the 
characteristic  function  for  a  normal  distribution,  using  the  trig  identity 
s;n2e  -  (!  -  cos26)/2,  and  again  using  Bessel's  integral  identity,  this 
time  for  Ed-  (15)  leads  to, 

E{  Jn(x)  }  =  exp[  -o2/4  ]  ■  ln/2(<72/4),  (25) 

' :r  n  even,  and  zero  when  n  is  odd.  This  expression  can  be  confirmed  by 
evaluating  the  expected  value  directly  with  the  help  of  an  integral 
tabulated  on  p.  710  of  reference  [6]. 

L  Let.  G(z)  be  the  generating  function  for  a  Poisson  distribution, 

G(z)  :  exp  [  rT  (z  -  i )  1. 


Putting  this  in  Eq.  (16)  and  integrating  yields, 

Ei  n!  /  sn+!  }  =  1  /(s  -  n  )  •  expl  -n  ].  (26) 

and,  in  particular,  when  s  =  1. 

E{  n!  }  =  l  /( 1  -  n  )  ■  exp  [  -n  ]  . 

This  is  easily  verified  to  be  correct,  as  well  as  the  fact  that  for  a 
Poisson  distribution  E{  n*  }  is  only  f mite  for  n  <  I . 

F\  !f  we  substitute  the  generating  function  for  a  Poisson  r.v.  into 
Eq.(18)  and  perform  the  integration,  we  easily  obtain  (letting  0=1), 

El  (2n  -  1)1!  }  =  I //(l  -  2nT  )  •  exp[  -n  ].  (27) 

Clearly,  this  is  finite  only  for  n  <1/2. 

Agam  using  the  generating  function  for  a  Poisson  r.v.,  Eq.  (20)  yields, 
after  using  a  trigonometric  identity  for  sin20  and  Bessel's  integral 
representation  for  the  Bessel  function  of  zeroth  order, 

El  (2m;n  /  (2m  *  1  )M  }  =  expl  -n/2]  J0(i  n/2) 

=  expl  -n  /  2  ]  l0(  n/2),  (28) 

where  In  is  the  modified  Bessel  function  of  zeroth  order. 


VjLL  Conclusion 

Some  of  the  foregoing  integral  identities  involving  character ist ic 
functions  and  generating  functions  may  be  derived  or  verified  using  other 
methods.  For  example,  I  originally  obtained  Eg.  (2)  by  expanding  Cx2(<3) 
in  a  McLaur in  series,  replacing  the  derivatives  wrt  %  by  even  order 
derivatives  of  Cy(0  wrt  £,  and  re-summing  the  infinite  series.  That 

derivation,  which  relies  on  the  existence  of  all  moments,  is  presented  in 
the  Appendix.  Similarly,  the  expression  for  E{  1  /(n  ♦  I)}  follows  easily 
from  integrating,  term  by  term,  the  infinite  series  definition  of  G(z).  In 
fact,  expressions  for  fractional  and/or  inverse  moments,  including  some 
of  those  derived  in  Section  IV,  have  been  expressed  elsewhere'’-8-9  in  a 
unified  manner  in  terms  of  fractional  mtegro-dif f erent lat ions  of  the 
MGF,  generalizing  the  usual  formulas  for  moments  and  factorial 
moments. 

However,  alternate  derivations  are  not  readily  identified  for  all  of  our 
■ntegral  relations.  The  point  is,  that  by  presenting  our  unified  treatment 
'.containing  as  a  proper  subset  some  of  the  previously  mentioned 
formalisms)  it  becomes  straightforward  to  obtain  new  integral 
dentines  for  random  variables  by  a  judicious  search  of  tables  of 
integrals  such  as  reference  [61.  As  a  final  example,  a  somewhat 
gratuitous  result  is  obtained  by  consideration  of  integral  no.  3  on  p.  304 
in  reference  [6], 

00 

j  exp  i  -pt  )  /[  I  ♦  exp(-qt)  1  dt  -  m/g  cosec  [ptr/q],  q  >  p  >  0, 

or  q  <  p  <  0. 


-00 


Let.  q  =  l,  p  ->  x.  a  r.v.  in  the  interval  [0,  I]  and  average  over  x, 

00 

1/tt  J  Nx(-t)  ^ 1  +  exp(-t)  ]  dt  =  E{  cosec  [tt-x]  }.  (29) 

-cc 

It  is  clear  that  many  other  integral  identities  for  random  variables  can 
be  generated  in  the  same  manner.  The  only  requirement  is  that  the 
implicit  interchange  of  orders  of  integration  be  justified. 


Acknowledgements:  The  author  is  grateful  for  advice  and  encouragement 
from  Profs.  E.  McKenzie,  D.  Gaver,  and  P.  A.  W.  Lewis;  and  to  Prof.  A. 
tawrance  for  references  [7,  8],  This  work  was  partially  supported  under 
j  contract  from  U.  S.  Navy,  SPAWARS,  PDW107-5. 


17 


REFERENCES 


' .  w.  Feller,  An  lot 


to  Pr  ( 


.ication; 


volume  II.  (Wiley,  New  York)  1965,  p.  9/7 


2.  w.  Feller,  An  Introduction  to  Probability  Theory  and  Its  Applications. 
Volume  I,  (Wiley,  New  York)  I960,  p.  298 

3.  E  Rockower  and  N.  Abraham,  "Calculating  Generating  Functions  from 
Characteristic  Functions,  with  application  to  quantum  optics", 

J.  Phys  A.  11,  i 0.  DD  1889  -  1879  (1978) 

9.  N.  A.  C.  Cress'e,  A.  5.  Davis.  J.  L.  Folks,  and  G.  E.  Policello  II,  "The 
foment -generating  Function  and  Negative  Integer  Moments", 

The  American  mat  sYcan  35.  No.  3,  pp.  148-150  (198!) 

5  M  mac  and  w  E  S’ rawderman.  "Negative  Moments  of  Positive 
Random-  var  ab.es".  .  Amer.can  Stat  Assoc  67,  No.  338, 


v  ;r  -  v  Rvm  *  'jbms  jf  integrals.  Series,  and 


'*  ‘  1  w,  T  ' 


<  c  ,r  •  ye,;. 

■;  M.  in  pni  preprint  C/ t»o.TS- 
-  '■  ^r"  -  -if  Dr-cif'w&  Qijgqrat'r 


No^a  v ar  aD;ec  appear  m  Austr  J  stat 


Y ,q  ^  5  „  pr" 


^ ament  Generating  cjnct:on  Has 


'.meats'.  2*  at  -st  F i arm  m-ference  13..  337  -  344  ( 19861 
vyj“-:r.  a  -'estise  or,  ~.ne  Tmeori  cf  Bessel  Functions.  2nd  ed 


Cambr-dqe  un ' ver s < f u  Press'  1 9S8.  p  389 

-1  AbramowfT  a^d  I  Stequn.  Handbook  of  Mathematical  fjjnctjons, 

■  r,r-  ,r  p  -,nc  ‘Jc,u,  ,  -,r  L  '  '  Ohc 


APPENDIX 


In  this  appendix  an  alternate  derivation  of  Eg.  (2)  is  presented,  patterned 


after  that  in  reference  [3]  in  which  it  is  applied  to  help  in  the 


interpretation  of  a  complicated  characteristic  function  and  identify 
independent  processes  which  combine  to  form  the  final  process.  The 
usual  definition  for  the  characteristic  function  yields, 


Cy(0  =  E{  exp(i£x) }  ,  and  Cx2(#)  =  E{  exp(i#x2)  }. 


We  solve  the  problem;  given  Cx  ,  to  find  Cx2  .  First,  note  that  the  even 


moments  of  x  can  be  expressed  alternatively  as, 


i  l/i2  (d/dO2  ln  Cx(0|  j=0  =  E(x2nl. 


!  i/i  (a/a?)  ]nC,.2(?)  „-n  =  E(  xzn  ) 


:.e.  the  !hs  of  these  two  expressions  are  equal,  assuming  the  moments 


e:cst. 


We  can  new  write  the  ordinary  McLaurin  series  for  Cx2  as. 


?('#)  --  2  [1/i  Cd/d^')  ]m  Cx2(^'){^-=0  /ml 


which  car.  re  rewritten  using  Eqs.  (A2)  as. 


cx 2(#)  =  2  [  -(d/dC)2  r  cx(0|  p=o  0?)m  /m! 


¥ 


Next,  re-sum  this  power  series  in  (d/d£)2  to  obtain, 


Cx2(2)  =  exp  {  -\Z  (d/d£)2  )  CX(C)|^=0  •  (A5) 

Now,  note  that, 

exp(  o2(d/dx)2 }  r(x)  =  1//(2ttc2)  j  exp{  -(x  -  x')2A2o2)  f(x’)  ax', 

which  can  be  verified  using  the  convolution  theorem  of  Fourier 
transforms.  This  equation  could  also  have  been  obtained  directly  from 
Eq.  ( 1 )  if  we  promote  x  to  be  the  operator  d/dx,  instead  of  a  r.v.,  and 
then  post-mu  1 1 ip iy  by  f(x).  Making  the  change  of  variables  x  -->  £, 
c2  -->  -2itf,  and  setting  £  =  0  we  recover  Eq.  (2). 
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